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Abstract 

We calculate the rate of the emission of the photons from the QGP. We base on 
the real-time kinetic approach [1] without an explicit assumption about a complete 
thermal equilibrium in the emitting system. 
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1. The subject of this note is as follows. In recent papers devoted to the rate of 
emission of the real photons from the quark-gluon plasma one may find very different 
final answers. Most of them agree as for the main logarithmic terms which appropri- 
ately estimate rate of emission only at unphysically large energies. Next terms which 
come to be significant at the typical value of E/T ~ 10 different authors write in 
manifestly different way [2,3]. It looks extremely surprising even the different func- 
tional dependence from the basic parameters, though all authors start from the same 
initial formulas. I have seen no papers were the procedure of analytic approximation 
would have been discussed in detail except paper [4] were expression for the rate of 
emission was given in the form of one-dimensional integrals which allow the reliable 
further analysis by very simple means. 

2. Expression for the rate of emission, more accurate then the main logarithm, is 
needed in connection with some modern scenario of the QGP formation which predict 
suppression of the quark component against background of the hot glue rather than 
complete thermalization of the QGP [5]. 

For sake of completeness we shall start with reminding of the definitions and the 
main framework of calculations. For photon spectrum the rate of emission is given 



In what follows we use the same way of derivation of the two loop expansion of 
7rfo as in [1,4] and denote g^ioi—k) as n(k). 

3. As the threshold of the first Born's term and the one-loop radiative corrections 
is k 2 > Am 2 they are absent for real photons. Physically, this mean that the two 
fermion states with the time-like momenta can not produce the massless photon. In 
the vertex which emits the photon one of the fermion legs should carry the time-like 
momenta and another - a space-like one. The fermion field excitations with space-like 
momentum can not belong to normal set of the normalizable states and therefore to 
contribute the density matrix. They can be only virtual. So only the real processes of 
the next perturbation order can contribute the rate of emission. The desired virtuality 
of the fermion line in the electromagnetic vertex comes from the intermediate state 
of the annihilation or the Compton process. 

It would have been inconsistent to include the Born's term with the fermion lines 
dressed by the one-loop self-energy in the whole scheme where the two-loop and the 
next graphs are considered perturbatively. It will lead to an over-count in the subset 
of graphs of the self-energy type and to the under-count in the group of graphs of the 
vertex type. 



by [1,4] 




(1) 



with the usual definition of two currents correlation function 




(2) 
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In fact, an account for some elements of a collective behavior in the ensemble of 
the partons is needed only in a very narrow region of a collinear geometry of the real 
processes which leads to the logarithmic divergence of the emission rate. A simple 
restriction of the size of the domain of the coherence of the interaction between the 
partons reflects the physics of the partonic ensemble and resolve the problem of a 
divergence. 

This allows us to consider only the real processes, annihilation and Compton, 
where the necessary virtuality is provided by the gluons. In this case the trace of the 
electromagnetic polarization tensor reads as [6] 

j p 2 „2 AT (~< r 

^eai = J d A pd A q5{q 2 )5[{p - k) 2 - m 2 }5[(p + qf]{SW a + SW C ) X 

x A h _g_ H C i 

(p 2 — m 2 ) 2 p 2 — m 2 p 2 — m 2 + 2kq 

where we denoted: 

A = 4m 2 , B = 3m 2 + 2(kq)-^--, C = m 2 + ^. (4) 

(kq) (kq) 

The expression in the curly brackets is a sum of the squared moduli of the matrix 
elements of the annihilation process, qq — > gj, or Compton process, qg — > qj and 
qg — > qj. Just for the purpose of the following calculations it is written not it terms 
of habitual Mandelstam variables (s, t, u). Specification of the process as well as their 
kinematic boundaries are due to statistical weights. 

In the formal approximation when distributions of initial particles are taken to be 
Boltzmann ones the statistical weight of the annihilation process looks as 

SW a = 9(k - p o )0(q o + Po)0(q o )n F (k o - p )n F (q + p )n B (q ) w 

e -{ku)/T 

w 9(k - p )9(q + p )9(q ) eHqu)/T _ x ■ (5) 

For the Compton rates the statistical weight is equal to 

SW C = -9(p + qo)9(p - ko)9(-q )n F (p + q )[l - n F (p - k )]n B (-q ) - 
-9(- Po - q )Q(k - p )9(-q )n F (k - Po)[l - n F (- Po - q )}n B (-q )- w 
» -e~ iku ^ T [9{p Q + q )9(p - k )9(-q )[ex V ((pu - ku)/T) + l]" 1 - 

-9(-po - qo)9(k -po)9(-q )[exp((qu)/T) + l]" 1 ] (6) 

This approximation was proved to be quantitatively effective for the case of the true 
thermal equilibrium of the emitting plasma. 
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We may also assume that quarks and gluons are not in thermal equilibrium. Com- 
puter simulation of the parton cascade usually produces something similar to the 
Boltzmann distributions damaged by the specific parameters (q and Cg- We adopt 
for quarks and gluons 

n F (p) = ( Q e- pu/T , n B (p) = C G e~^ T , (7) 

where u^, the 4-velocity of the nonequilibrium partonic media, fugacities ( and tem- 
perature T are very smooth functions of space-time coordinates. Scenario of the hot 
glue leads to (q < (a < 1 an d all these quantities have no habitual thermodynamical 
meaning. 

In this parametrization of the partons distributions which we shall name the par- 
tons cascade, the statistical weight of the annihilation process with emission of a 
gluon looks as 

sw casc =w 6{kQ _ po)e{go + po )0 ((?o ) C 2 e -(M/T (e -M/T + ^^W/T) (g) 

For the Compton rate of the photon emission the statistical weight equals to 

SWr° » -CqCg{0( Po + q )9( P o - k )6(-q )e-^ T [l - ^-^-ku) /T] + 

+0(- Po - qo)0(k o - Po )9(-q )e~ ku / T e^ u+pu ^ T [l - C Q e {qu+pu)/T }} (9) 

For any kind of these distributions we can perform an exact integration over p using 
the Breit reference system where k + q = 0, c.m.s. of the reaction qq — » g^y. 
For the equilibrium distributions of the partons this integration leads to 



ie 2 - 2 



9 -^e-^\ d* q 5( q y( qo )e[(k + q ) 2 -4 m 2 }£^ (10) 



where 



„ 2m 2 2m\, l-Jl-2m 2 x 2m 2 , 2m 2 
F a {x) = 1 + In V + 1 + Wl 11 

x x 2 l + ^l-2m 2 /x % V x 

For the Compton process it is reasonable to start with a chain of changes of 
variables, p — > — p + k — q, in the first term and q — > — q — p in both terms. After 
that we may also easily perform the integration over p using the same Breit reference 
system. It gives 



T^compt 



.*£^^jW-^)J*!!L (12 ) 
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where q is the momentum of the (anti) quark in the final state and 

^, N 2m 2 2m 4 , , m 2 + 2x 4m 2 2x(m 2 + x) /1oX 
^ c x = l r In — + + - 1 J 13 

In course of these calculations we could easily learn that the energy p of a virtual 
quark in the annihilation process equals to zero in the Breit system. In the Compton 
process it equals to the ratio m 2 /2(ko + go) which is small for the hard photon also. 
At least in both cases the momenta of virtual quark are space-like. Because of the 
high energy of the photon the 4-momentum p^ of a virtual quark is very close to the 
light cone in the rest frame of the media. 

4. Further integration is somewhat tricky. Indeed, for the photons n compt depends 
only on one time-like 4- vector u, u 2 = 1. So to continue calculations in a reasonable 
way we are to chose reference frame u — 0. In this frame the integral of Eq.(9) takes 
form, 

- nT J™ d[ln(l + e-V " 2+m2 )]g J ^dzF c (k^q 2 + m 2 + \k\qz) (14) 

where ko — |k| stands for ku. The trick which allows one to reduce this double integral 
to a simple quadrature is as follows [4]: 

i) change the variable z for the new one, \k\qz = 77; 

ii) integration with respect to q by parts; 

iii) successive changes of variables: q = msinh^, then me^/T = 2y. 
The result reads as 

7W = >y e -^ J™ dy ln(l + e^-i) x (15) 

x[{1 -y-W )Hl + T ) + J + ltTW ] 

where £ = m 2 /{ku)T , A = m/T and singular point y = corresponds to the infinite 
rapidity (f) of the massless quark. 

For the qq annihilation into gluon and real photon, k 2 = 0, the angular integration 
is also performed in the rest frame of the media and may be reduced to 

ttT j d[\n(l - e- q/T )]q J ^ dz9(q(k + \k\z) - 2m 2 )F a ( y (k + \k\z)q) (16) 

After changing variable z for r\ = q(k + \k\z) integration over q by parts and changing 
q for y = (m 2 jku)q we arrive at 

ie 2 g 2 N c C Fr 



2vr 3 

y V 



roc 

_ T 2 e -ku/T J dy\n(l-e- y )x 




x[(! + 7. ~ h^ cosh ~\lT + (1 + 7.) \/l " '1 (17) 
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It is easy to see that the spectrum of gluons in the annihilation process is effectively 
cut from the above at gluon momenta ~ T. So from a simple kinematics it follows 
that at high photon energy, E » T, the energies of the initial quarks should be 
large, time-like and slightly above the light cone. 

The Compton process can be more precisely identified as the bremsstrahlung of 
a photon due to the hard scattering of a quark on a gluon. Maximum rate of the 
emission is produced by the infinite rapidity of the massless quark in the final state. 

5. Both annihilation and Compton rates as they are given above logarithmically 
diverge when quark mass is going to zero and the geometry of collision is collinear. 
In this case all the participants of the reaction interact infinitely long. This can not 
be true in the media and all the momenta should be cut off from the below. 

The external cut-off should come from the smallest of the three lengths. The first 
one is the Compton wave-length l c ~ 1/m. The second one is the mean free path 
defined via the ordinary or transport cross-section, l m f p ~ l/g 2 T (or ~ l/g 4 T). The 
third one is connected with the amplitude of the forward scattering, lf s ~ 1/gT ~ 
l/m>therm- As long as we keep quarks massless and consider coupling as the small 
parameter we must chose the length lf s as the smallest one. So whenever we meet 
the collinear singularity the momenta of both real and virtual quarks should be cut 
off from the below at ~ gT. 

For the production of photons with the energy E much exceeding any of these 
scales the particular choice is not so significant within the reasonable accuracy. The 
cut-off mass appears only under the logarithm and its variation is not so valuable. 

6. Having reached our goal to present the expression for the rates of emission 
from the equilibrium plasma in the form of simple quadratures we can easily find 
them numerically. We can analyze the reliability of different analytic approximations 
and compare them with those known in the literature. 

The way of calculations chosen in the papers [2] and [3] had followed paper [7] 
where at the very initial stage it was performed change of variables from original mo- 
menta to the invariant variables (s,t,u). This trick produces a set of delta- functions 
which are naturally used to get rid of angular integrals. This results in a cumbersome 
four-dimensional integral which can be simplified only by assuming that s » t. This 
assumption (as it was used in [7]) do not touch the main logarithmic term which is 
contributed basically by the low-angle scattering (small t). But it violates the con- 
stant term which is mainly due to large-angle scattering, not to tell about the next 
terms. Though these are parametrically small, they are not small numerically in the 
physically actual region. 

I failed to trace out the way of analytic approximation in Ref. [3] but the final result 
disagree even with the well established functional form of the leading logarithmic 
terms. 

Let us start with the annihilation mechanism of the real photons production. 
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Notice that rate of photon emission depend upon a single parameter £ = m 2 /{ku)T 
and contains the logarithmic (collinear) divergence when m — > 0. Assuming the 
regularization by means of the mass cut-off at m ~ gT we see that £ << 1 at least 
because of (ku) » m ~ T. 

My statement is that the approximation which picks up only the main logarithm of 
the small parameter £ fails at reasonable values of this parameter (actually, £ ~ 0.1). 
I will show that terms like £ln£, £ln 2 £,... appear as well and that at reasonable 
values of £ they essentially change the rate of emission. 

Let us consider the leading integral of n anni 

- dy\n{l - e~ y )\n^- (18) 

Its main part which is easily obtained by taking zero for the low limit equals to 

7r 2 4 , r°° y In y dy . . 

ir< ln i + 1) + { /» W K -°- 24 > (19 > 

while estimation of the residue at small £ adds a term like 

^ln^-£) (20) 

Estimation of the next terms is somewhat more complicated but it is very easy to 
understand the type of corrections: as the integrand has logarithm in the nominator 
so the extra powers of y in the denominator will produce the powers of ln£. 

Just similar analysis can be performed for the Compton rate of emission. The 
integral (15) is effectively cut from the below by y min ~ A 2 /4 . The main term may 
be approximated as 

jTM^ + l)0.| + i) (21) 
Proceeding as above we get the leading part, 

^(ln | + I) + { jH ln(l + e~y) \nydy = -0.372} (22) 

and the corrections, 

_^(l n ^_l)l n2 + !L^i n 2 (23) 

4 v 4 ; 12 8 K J 

The "reasonably small" value of the parameter A = m/T is 0(1) and the correc- 
tions are comparable with the main term at the previous value of £ ~ 0.1. The rest 
terms in the integral (15) will produce corrections with the next powers of In A, just 
as in the case of the annihilation . 
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This kind of analysis helps us to understand why it may be dangerous to use the 
limit of the high photon energy in the many- dimensional integrals. In this form it is 
hard to figure out the true structure of corrections. 

At this circumstances it seems more reliable to use the quadratures (15) and (17) 
as the definitions without further approximations. The result of numerical calcula- 
tions is plotted at Fig.l apart from the common constant factor and the Boltzmann 
factor exp(—E/T). The difference between exact function and the approximate ones 
is easily seen to be almost constant at the photon energies exceeding 2 Gev. An 
approximation given in Ref.[2] approaches the exact result only at E ~100Gev. At 
E <4Gev this approximation leads even to the wrong relation between the annihila- 
tion and Compton rates. 

7. The partons distributions (7) of the nonequilibrium regime lead to another ex- 
pressions for the rate of emission. All technique of calculations remains the same and 
the final results immediately follow from the next very simple observation. Inspecting 
statistical weights given by Eqs. (5) and (8) we see that the following changes are 
required in the annihilation channel: 

[e WT e -Wr ^ C 2 Q [^ qu/T + Cce~ 2qu/T ] (24) 

The middle of this string corresponds to the straightforward Boltzmann approxima- 
tion which results in the replacement of the ln(l — e~ y ) in Eq.(17) by the — e~ y . 
Now the prescription is evident, 



7T = 

ann 



ie g N c C F T2e _ ku/T ^ joo ^ + C _^ e _ 2y] 



2tt 3 " ^ Jt 

y %y 




x[(i + i - ^)2cosh- 1 J y - + (i + i) - i] (25) 



The first exponent in the integrand corresponds to spontaneous emission of a gluon, 
the second exponent with the extra factor (a corresponds to the induced emission. 

The rate of photons production in the Compton channel modifies in a similar way. 
A string of changes in the Eq.(12), 



results in 



„casc 
compt 



[ e WT + r i _> e -WT _> CqCg [ c -Wt _ (oe -¥/T] (26) 



x[(l-i-i; )ln(l + * + ^ + Mi±M] (27) 

y W i' y (e + %) 2j 
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where the second exponent is due to the Pauli suppression of a quark in the final 
state. 

The actual values of the fugacities which are the measure of the chemical equilib- 
rium are about (q ~ 0.5 and ( G ~ 0.75 [5]. In order to find out the influence of the 
kinetic nonequilibrium produced by the Boltzmann distributions we plot the rates 
given by (25) and (27) using (q = (c = 1 ( also apart from the common Boltzmann 
factor exp(-E/T)). 

We consciously do not integrate the rates over the history of the system which can 
be easily done for a simple form of the hydrodynamic background. The very approach 
of nonequilibrium field kinetics [1,6] was designed in order to use the nonequilibrium 
partons distributions generated in a "realistic" cascade. 
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